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ABSTRACT 

The apparent spectral evolution observed in the steep decay phase of many GRB early afterglows 
raises a great concern of the high-latitude "curvature effect" interpretation of this phase. However, 
previous curvature effect models only invoked a simple power law spectrum upon the cessation of the 
prompt internal emission. We investigate a model that invokes the "curvature effect" of a more general 
non-power-law spectrum and test this model with the Swift/XRT data of some GRBs. We show that 
one can reproduce both the observed lightcurve and the apparent spectral evolution of several GRBs 
using a model invoking a power-law spectrum with an exponential cut off. GRB 050814 is presented 
as an example. 

Subject headings: gamma-rays: bursts 



1. INTRODUCTION 

Most of the early X-Ray afterglows detected by Swift 
(Gehrels et al. 2004) show a steep decay phase around 
100^1000 seconds after the burst trigger (Tagliaferri et 
al. 2005). The main characteristics of this steep decay 
phase include the following. (1) It connects smoothly to 
the prompt 7-ray light curve extrapolated to the X-ray 
band, suggesting that it is the "tail" of the prompt emis- 
sion (Barthelmy et al. 2005, O'Brien et 2006, Liang et 
al 2006). (2) The decay slope is typically 3 ~ 5 when 
choosing the GRB trigger time as the zero time point 
t Q (Tagliaferri et al. 2005; Nousek et al. 2006; Zhang 
et al 2006). (3) The time-averaged spectral index of the 
steep decay phase is much different from that of the later 
shallow decay phase, indicating that it is a distinct new 
component that is unrelated to the conventional after- 
glow components (Zhang et al 2006; Liang et al. 2007). 
(4) Strong spectral evolution exists in about one third of 
the bursts that have a steep decay phase (Zhang et al. 
2007, hereafter ZLZ07; Butler & Kocevski 2007; Star- 
ling et al. 2008). All these features suggest that the 
steep decay phase holds the key to understand the con- 
nection between the prompt emission (internal) phase 
and the traditional afterglow (external) phase. Any pro- 
posed model (see Meszaros 2006; Zhang 2007 for reviews) 
should be able to explain these features. 

The so called "curvature effect", which accounts for 
the delayed photon emission from high latitudes with re- 
spect to the line of sight upon the abrupt cessation of 
emission in the prompt emission region (Fenimore et al. 
1996; Kumar & Panaitescu 2000; Dermer 2004; Dyks et 
al. 2005; Qin 2008a), has been suggested to play an im- 
portant role in shaping the sharp flux decline in GRB 
tails (Zhang et al. 2006; Liang et al. 2006; Wu et al. 
2006; Yamazaki et al. 2006). In the simplest model, it is 
assumed that the instantaneous spectrum at the end of 
the prompt emission is a simple power law with a spec- 
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tral index f3. The predicted temporal decay index of the 
emission is (with the convention F v oc t~ a v~P) 

a = 2 + /3, (1) 

if the time origin to define the log — log light curve, to, is 
taken as the beginning of the last emission episode before 
the cessation of emission. Adopting a time-averaged /? 
in the tails, Liang et al. (2006) found that Eq.fl]) is gen- 
erally valid. The strong spectral evolution identified in a 
group of GRB tails (ZLZ07) apparently violates Eq.([T]), 
which is valid only for a constant p. ZLZ07 then investi- 
gated a curvature effect model by assuming a structured 
jet with varying f3 at different latitudes and that the line 
of sight is near the jet axis 4 . One would then expect 
that Eq.fl]) is roughly satisfied, with both a and (3 being 
time-dependent. ZLZ07 found that this model does not 
fit the data well. 

These facts do not rule out the curvature effect inter- 
pretation of GRB tails, however. This is because the 
instantaneous spectrum upon the cessation of prompt 
emission may not be a simple power law. If the spec- 
trum has a curvature, as the emission from progressively 
higher latitudes reach the observer, the XRT band is 
sampling different segments of the intrinsic curved spec- 
trum (Fig.l). This would introduce an apparent spectral 
evolution in the decaying tail. The main goal of this pa- 
per is to test this more general curvature effect model 
using the available Swift XRT data. 

2. CURVATURE EFFECT OF A NON-POWERLAW 
SPECTRUM 

We consider a general non-power-law spectrum in the 
form of 

F„(u) = F v>c G(v) , (2) 

where G{y) is the function form of the spectrum with a 
characteristic frequency v c so that G(v c )=l, and F v>c = 
F u (v c ) is the normalization of the spectrum at v = v c . 

The curvature effect states that given a same spectrum 
at different latitudes with respect to the line of sight, 

4 Notice that this structured jet model is different from the tradi- 
tional one that invokes an angle-dependent energy/Lorentz factor, 
but not the spectral index (Zhang & Meszaros 2002; Rossi et al. 
2002). 
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TABLE 1 

Best-fitting parameters and their 1-sigma errors for the 
cutoff power curvature effect model for. grb050814. 



*0< s ) 



0.67(0.12) 10.2(1.3) 



1.56(0.25) 103.5(3.4) 0.002(0.04) 1 (fixed) 



one has F v ^ c oc V 2 and v c oc V, where T> is the Doppler 
factor. If the high-latitude angle 6 ^> T, the Dopper 
factor P oc T 1 , so that F v ^ c oc t~ 2 , v c oc t^ 1 (Kumar & 
Panaitescu 2000). Considering the tp effect (Zhang et al. 
2006; Liang et al. 2006), this can be written as 



and 



F„, c (t) = F v 



Vc{t) 



t-tp 

t„ — tp 



t-tp 

tp — tp 



(3) 



(4) 



for t 3> t p , where to refers to the time origin of the last 
pulse in the prompt emission and t p is the epoch when the 
curvature-effect decay starts (or the "peak" time of the 



lightcurve), F l 



F v ,c(tp) and v < 



c, P - v c {t p ). Notice 
P (a pure power law 
P- 2 . This is the 



that in the case of G{v) = (y/v c ) 
spectrum), one derives F v oc (t — to) 
relation Eq.fl]). 

We consider several physically motivated non- 
powerlaw spectra with a characteristic frequency v ci in- 
cluding the cut-off power law spectrum and the Band- 
function (Band et al. 1993). To explore the compati- 
bility with the data, we also investigate different forms 
of the cutoffs with varying sharpness. In all cases, the 
F Vp (t) and v p (t) follow Eqs.© and ®. When v c {t) drops 
across an observational narrow energy band, e.g. the 
Swift /XRT band, it introduces an apparent spectral soft- 
ening with time, which, if fitted by a power law, shows an 
increase of photon index with time. In the meantime, the 
flux within the observing band drops down rapidly, lead- 
ing to an apparent steep decay phase in the lightcurve 
(Fig.l). 

3. DATA REDUCTION AND SIMULATION METHOD 

We consider a time-dependent cutoff power law photon 
spectrum taking the form of 



N(E,t)=N (t) 



— Y 

1 keV/ 



exp 



E 



E c (t) 



(5) 



where r = (3 + 1 is the power law photon index, and k 
is a parameter to define the sharpness of the high energy 
cutoff in the spectrum, E c {t) is the time-dependent char- 
acteristic photon energy, and AT (t) is a time-dependent 
photon flux (in units of photons ■ keV -1 cm~ 2 s -1 ) at 1 
keV (Arnaud 1996). The choice of this function was 
encouraged by the fact that the spectral evolution of 
some GRB tails can be fitted by such an empirical model 
(Campana et al. 2006; ZLZ07; Yonetoku et al. 2008). 
According to Eqs.Q and (|4]), and noticing the conversion 
between the photon flux and the emission flux density, 
i.e. F v oc EN(E), we get 



E c (t) = E c 



t-tp 
t p — <q 



(6) 




Fig. 1. — A schematic picture showing that shifting a set of 
non-power-law spectra in time can equivalently give an apparent 
spectral evolution in a fixed band. The dashed lines represent a set 
of exponential-like spectra, whose F u (t) and fp(i) drop down with 
time according to Eqs. ((3j and (|4jl. The two vertical solid lines 
bracket the observed energy band. The thick solid lines denote the 
effective power law fits to the time-dependent spectra at each time 
step. 
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Fig. 2. — The lightcurve (upper panel) and spectral evolution 
(lower panel) of the X-ray tail of GRB 050814 with the best-fit 
theoretical model (black curve in upper panel and green curve in 
lower panel). The blue and red data points are the window timing 
and photon counting data, respectively. The inset shows time- 
dependent theoretical spectra with the XRT band (0.3-10 keV) 
bracketed by two vertical lines. The integers denote the time seg- 
ments for the time-resolved spectral analysis. 



and N{E C , t) = N c , p [(t - t )/(t p - t )] \ where A c , p 
N(E c ,t p ), and E CiV = E c (t p ). This gives 



iV (f) = N , 



t-tp 
t p — <o 



-U+r) 



(7) 
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Fig. 3. — Upper panel: Examples of simulated time-dependent 
spectra of GRB050814 with the best-fit parameters. The time in- 
tervals are 1,6 repectively as denoted in Fig. 2. In each panel, the 
data histogram displays the simulated spectrum, and the solid line 
displays the best-fit (\ 2 /dof = 39.0/61, 25.2/25) power law model 
(wabs * zwasb * powerlaw in XSPEC) that is used to derive the 
time-dependent photon index T. Lower pabel: The corresponding 
observed spectrum in the three time intervals and their power law 
fits (x 2 /dof = 47.1/46, 22.0/19). 
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Fig. 4. — The simulated cut-off power law spectrum at t p = 144 
s based on the best fit model confronted by the BAT data in the 
time interval (141.5 — 146.5) s. The reduced \ 2 of the fitting is 1.2 
with dof= 197. 

Notice that t p is the beginning of the steep decay, which 
is a parameter that can be directly constrained by the 
data. For a complete lightcurve, we read t p off from the 
lightcurve. In the case of an observational gap, usually 
t p can be reasonably fixed to the end of the prompt emis- 
sion. We therefore do not include this parameter into the 
fits, and derive the other five parameters, namely, No tP , 
E CjP , r, tg, and k from the data. At any time i, the model 
spectrum can be determined once these parameters are 
given. One can then confront the model with the real 
GRB data. 

The procedure includes the following steps. (1) For a 
given burst, we extract its Swift/XRT light curve and 
n slices of time-dependent spectra using the standard 
HEASoft/Swift Package. The details of the data reduc- 
tion method were described in ZLZ07. (2) Given a trial 
set of parameters in the theoretical spectra 5 {No,p, E CtP , 
F, to}, using Eqs.([5][5]) we model n time-dependent the- 
oretical spectra that correspond to the time bins that 
are used to derive the time-dependent observed spec- 

5 Notice that k is fixed to a certain value for a particular model, 
and is varied when different models are explored. 



tra. (3) Based on the theoretical spectra of each time 
slice, we simulate the corresponding model spectra by 
taking account of the observational effects, including the 
Swift /XRT response matrix, the absorption column den- 
sities (Nji) of both the Milky Way (extracted from the 
observations from step 1) and the host galaxy of the burst 
(a free parameter), the redshift (if applicable), and a 
Poisson noise background. Notice that riH,host is another 
parameter introduced in the model spectra (besides the 
other parameters introduced in the theoretical spectra). 
All these faked spectra can be obtained using HEASoft 
(Version 6.4) and Xspec (Version 12.4) (4) We fit the 
faked model spectra with a simple power law model, i.e. 
wabs * wabs * powerlaw (or wabs * zwabs * powerlaw if 
the redsift is available) in XSPEC and get the simulated 
fluxes and spectral indices of the n slices. Here the col- 
umn densities of both the Milk Way and the host galaxy 
are fixed to the obsered values as in Step 1. (5) We com- 
pare the simulated fluxes and spectral indices with the 
observed ones and access the goodness of the fits using \ 2 
statistics. (6) We refine the trial set of parameters based 
on the comparison and repeat steps (2)- (5) when neces- 
sary. We test whether we can reach a set of best-fitting 
parameters that can reproduce both the light curve and 
the apparent spectral evolution as observed. 

4. AN EXAMPLE: GRB050814 

We apply the method to GRB050814, a typical burst 
with well-observed X-ray tail with strong spectral evolu- 
tion. As seen in Fig. 2, the tail has a steep decay index 
of ~ 3.2, and a strong spectral evolution is apparent at 6 
t < 600 s. These features are common in most of the 
GRB X-ray tails. We first fix k = 1 in Eq. [5j which 
corresponds to the simplest cutoff powerlaw model. The 
initial trial parameters we choose are (r, No jP , to, E Pj q, 
n H ,host) = (1.2,0.4,72.0,30.0,0.05). The peak time t p 
is fixed to 143.6 s, which corresponds to the end of the 
prompt emission. Some IDL scripts are developed to fol- 
low the procedure described in §3 to automatically search 
for the best-fit parameters to match both the observed 
light curve and the time-dependent spectral index. The 
final best-fitting parameters are shown in Table 1. The 
corresponding simulated light curve (black curve) and 
spectral indices (green curve) are shown in Fig. 2. Figure 
2 suggests that the sharp decay and the spectral evolu- 
tion in the tail of GRB 050814 can be indeed explained 
by the curvature effect with a cutoff power law spec- 
trum. In Fig. 3 we present the comparison between the 
simualted and observed spectra in the time steps 1 and 
6 (as examples), which show reasonable consistency. 

Our model predicts that the prompt emission spec- 
trum at t p ~ 144 s should be a cut-off power law with 
the parameters in Table 1. In order to confirm this, we 
subtract the BAT-band spectrum in the time interval 
(141.5 — 146.5) s, and compare with the data with the 
model prediction. As shown in Fig. 4, the BAT data is 
roughly consistent with the model prediction, suggesting 
the validity of the model. 

Some physical parameters can be constrained accord- 
ing to our model. The time interval from t p to the begin- 

6 The PC mode spectra become harder at the end of this tail, 
which might be due to the contamination of the harder shallow 
decay component. For simplification, we focus on the WT mode 
data only. 
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ning of the steep decay phase tt a a,o may be related to the 
angular spreading time scale T ang — (t ta ii.o - t p )/ (1 + z). 
Noticing z ~ 5.3 for GRB050814 (Jakobsson et al. 2005), 
we can estimate the Lorentz factor of the fireball as T — 
(R/2cT ang ) 1 / 2 ~ 69i?^ 2 , where R 15 = i?/(10 15 cm) is 
the normalized emission radius. Since we know the spec- 
tral peak energy E p at t p , we can also estimate the corre- 
sponding electrons' Lorentz factor for synchrotron emis- 

sion by 7e , p = [E p /(hT^ c )] 1/2 ~ 2.4 x lO 3 ^ 4 ^ 1 / 2 . 
From the rest frame duration of the X-ray tail we are an- 
alyzing Ttail = {ttaU,e-ttail,0)/0-+z) ~ (378-165)/6.3 = 

33.8 s, one can constrain the minimum jet opening an- 
gle as 9j > {2cT tai i/R) 1 / 2 = 2.6° x R^ /2 . These values 
are generally consistent with those derived from various 
other methods. 

We find that the abruptness parameter k cannot be 
very different from unity. A Band-function spectrum in- 
troduces a less significant spectral evolution and it can- 
not reproduce the data (cf. Qin 2008b). 

5. DISCUSSIONS AND CONCLUSIONS 

We have successfully modeled the lightcurve and spec- 
tral evolution of the X-ray tail of GRB050814 using the 
curvature effect model of a cutoff power law spectrum 
with an exponential cutoff (k = 1). It has been discussed 
in the literature (e.g. Fan & Wei 2005; Barniol-Duran & 
Kumar 2008) that the GRB central engine may not die 
abruptly, and that the observed X-ray tails may reflect 
the dying history of the central engine. If this is in- 
deed the case, the strong spectral evolution in the X-ray 
tails would demand a time-dependent particle accelera- 



tion mechanism that gives a progressively soft particle 
spectrum. Such a behavior has not been predicted by 
particle acceleration theories. Our results suggest that 
at least for some tails, the spectral evolution is simply a 
consequence of the curvature effect: the observer views 
emission from the progressively higher latitudes from the 
line of sight, so that the XRT band is sampling the dif- 
ferent segments of a curved spectrum. This is a simpler 
interpretation. 

The phenomenology of the X-ray tails are different 
from case to case (ZLZ07). We have applied our model 
to some other clean X-Ray tails, such as GRB050724, 
GRB080523, and find that they can be also interpreted 
by this model. Some other tails have superposed X-ray 
flares, making a robust test of the model difficult. A 
systematic survey of all the data sample is needed to 
address what fraction of the bursts can be interpreted in 
this way or they demand other physically distinct models 
(e.g. Barniol-Duran & Kumar 2008; Dado et al. 2008). 
This is beyond the scope of this Letter. 
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